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FSI:  STATE  OF  THE  ART 


Structural  dynamics 


-  multibody  dynamics 

-  geometrical  nonlinearities  (large  displacements,  rotations 

and  strains) 

-  material  nonlinearities  (nonlinear  constitutive  models) 

-  crack  propagation  (failure) 


MIIlM 


nwt.i 


-  shocks 

-  turbulence 


-  static  (steady)  and  dynamic  (unsteady) 


-  eigen 


IMPACT  ON  ENGINEERING 


Strong 


-  analysis  of  carefully  selected  critical  configurations 


Weak  (or  not  as  strong) 


-  routine  analysis 

-  design  (and  test)  [>-  significant  CPU  time  issues 


-  control 


“  [If  I  am  not  getting  the  NASTRAN 
answer  after  4  hours  on  a  Cray, 
then  God  is  sending  me  the 
message  /  have  the  wrong  design ]  " 

Burt  Rutan,  1993 


Manuel  Balce  Ceneta/  AP 


REDUCED  ORDER  MODEL  (ROM) 


Model  reduction  (MOR) 


-  build  the  lowest  dimensional  model  that  can  capture  the 
dominant  behavior  of  the  system  of  interest  by  projecting 
the  high-fidelity  model  onto  a  well-chosen  subspace 

i=>  drastic  CPU  time  reduction 


Complex,  time-dependent  problems  (with  a  CFD  component) 
o  Perturbation  problems  (stability,  trends,  control,  etc.) 


-  linearized  1 - !>  linear  ROMs 


o  Response  problems  (behavior,  performance,  etc.) 


-  nonlinear  1 - !>  linearized  1 - v>  linear  ROMs 

Newton 


MODEL  REDUCTION 


What  this  is  NOT  about 


-  building  the  simplest  model 

-  building  a  variable  (or  multi)  fidelity  model 

-  adopting  the  coarsest  mesh 

-  substructuring 

-  constructing  a  "surrogate"  or  "meta"  model 


LINEARIZED  DYNAMICAL  SYSTEMS 


input  output 

- *  I  SYSTEM  - ; 

M- —  [Mi#  *  •  •/M'p]  y  —  [yi**  •  'lYp  ] 


External  description  (input-output  map) 


y  =  S(u)  =  h*u  =  L^h(t-  t)  h(t)  dx 


Internal  description  (model) 


u  =  A(m,  •••<  MP)u  +  f 

y  ^ ( u ,  t,  m,  •  •  f-ip) 

u(t0)  =  u0 


SURROGATE  MODEL 


I 


External  description 


input 

[Hi.---.mJ 


fl 


SURROGATE  °utput 

MODEL  y  =  [yi.  -  ,yp] 


High-fidelity 

analysis 

Surrogate 

analysis 


Input 


-  meta-model,  surrogate  model,  response  surface,  kriging 


-  lower-dimensionality  is  not  guaranteed 

-  it  is  not  a  model  of  the  system  but  of  the  output 


MODEL  REDUCTION 


Internal  description 


=  A 


Mi»  •  •  •'  M 


)u  +  f  (dimension  =  N) 


rojection  onto  a  subspace  of  dimension  k  <<  N 


-  right  Reduced-Order  Basis  (ROB)  ,  k  <<  N 


u  -  'ky 


Vky  =  A( 


)  ky  +  f  +  r 


-  left  Reduced-Order  Basis  (ROB)  ,  k  <<  N 


constraints 


y  -  utTA( 


)Vty  +  lVf 


parameterized  linear(ized)  Reduced-Order  Model  (ROM) 


MODEL  REDUCTION 


issues 


choice  of  the  lower  dimension  and  the  ROBs  and 
dependence  of  the  resulting  ROM  on  the  parameters  { 


is  a 


LINEAR(IZED)  FLUID-STRUCTURE 


General  purpose  linearized  aeroelastic  high-fidelity  model  (HFM) 


MULTIDISCIPLINARY  ROM 


CFD-based  HFM 
Balanced  POD-based  ROBs 
(with  stability  guarantee) 


FE-based  HFM 
Eigen-based  ROB 
(with  truncation) 


1 


FLUID-STRUCTURE  ROM 


Fluid  Structure  Coupling 


ROBs:  Uk,  Vk  /  UkTVk  =  I  ROB:  Xm  (XmTMXm  =  l)  Br=UkTBXm 

w  =  Vkwr  u  =  Xm  um  Cr=  UkT  C  Xm 

ROM:  Hr  =  UkT  H  Vk  ROM:  Qm2  =  X^  K  Xm  Pr  =  Xj  P  Vk 


General  purpose  aeroelastic  ROM  (flutter,  response,  control,  ...) 


OFFLINE  FLUID  SNAPSHOTS 


Natural  mode  shapes  (ground  vibrations) 


9  modes) 


Excitations  in  a  sampled  frequency  range 


0  <  K 


nf  =  5  frequencies 


Responses  of  linearized  flow  (frequency  domain) 


H  +  ia>J 


(C  +  i©„B)  uj  j  =  1 


data  compression 

primal  fluid  snapshots  W  - * 


(-HT  +  i©  I) 


data  compression 

dual  fluid  snapshots  W* - *• 


DATA  COMPRESSION 


Modal  superposition  (Fourier  decomposition) 


-  limited  range  of  applications 


Proper  Orthogonal  Decomposition  (POD) 


-  lacks  stability 


Balanced  Proper  Orthogonal  Decomposition  (POD) 


-  more  robust  than  POD  but  still  lacks  stability 


GUARANTEED  STABILIZATION  METHOD 


ROM  stabilization  method  ( 


Amsallem  and  Farhat,  2011 


-  universal 

-  non  intrusive 

-  computational  complexity  scales  with  the  size  of  the  ROM 

-  preserves  accuracy  of  original  (unstable)  ROM 


input:  ROM  Stabilization  outPut:  stable  R0M 

Method 


F-16  BLOCK  40 


Higher-fidelity  (higher-dimensional)  models  (HFM) 

-  structure  :  FEM  with  168,799  dots 

-  fluid  :  Euler  CFD  model  with  403,919  grid  points 

ROMs 

-  structure  :  projection  on  ROB  consisting  of  first 

9  natural  mode  shapes 

-  fluid  :  projection  on  ROB  of  dimension  60  generated 

by  POD  using  99  snapshots 

(9  shapes  x  5  cos  x  2  +  9  shapes  at  0  Hz) 


Lift  (Ibf) 


PARAMETER  SENSITIVITY 


POD-based  fluid  ROB  (60)  built  at  M  =  0.710  (trimmed  angle) 


L-HFM :  M=0.710 
(2,188,394  dofs) 

ROM  :  M=0.710 
(69  dofs) 

L-HFM :  M=0.799 
ROM  :  M=0.799 


Time  (sec) 


ANALOGY 


Configuration  1  Configuration  2 


Reduced-order  basis  Vkl 


Reduced-order  basis  Vk2  =  Vkl 


TURNAROUND  TIME 


F-16  Block  40  —  1  operating  point  —  2nd-order  discretization 


64-processor 


Linux  cluster 


Construction  and  processing  of  aeroelastic  ROM  in  TO,  1.0 1  s 


rnmm j-f ^Tmn 

I  tion  of  99  fluid  snapshots  i 

iction  of  fluid  ROM 

Processing  aeroel&st+c-ROM 

3U  time 


50  minutes 

0.25  minu 

0.10  minute 
57.25  minu 


Parameter  2 


APPROACH 


Database  of  fixed-size  stable  ROMs 
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□  : 


Precomputed 


ROM 

Instantaneous 

ROM 

Zonal 


interpolation 


Parameter  1 


CONSISTENCY 


0i - 1 - 1 - 1 - 1 - 1— 

0  C.5  1  1.5  2  2.5 


Eigenvectors  4-  -j-  X  ~  ~f"  ~f ''  “f-  v2 

precomputed  ROMs  are  not  necessarily  computed 
in  a  consistent  set  of  generalized  coordinates 


ROM  TRANSFORMATION 


Congruence  transformations 


-  reference  ROB:  Vref  =  Vk(juiref) 

-Vi  * - VjQi  Qi  =  arg  min  ||  V;  Q.- Vref||F 

Q.eO(k) 


Orthogonal  Procrustes  problem 


-  SVD:  ViTVref  =  U.S;  Z,T 


Qi=  UjZj 


QUOTIENT  AND  EMBEDDED  MANIFOLDS 


-  standard  interpolation 
of  a  ROM  does  not 
produce  a  ROM 

-  interpolation  must  be 
carried  out  on  a  manifold 
characterizing  the  ROM 


Manifolds  of  interest  (quotient  =  blue,  embedded  =  yellow) 


-  span(Vk)  belongs  to  the  Grassmann  manifold  G(k,N) 

-  Ar(  )  =  VkTA(  )Vk  belongs  to 

*  manifold  of  invertible  matrices  GL(k)  [fluid] 

*  manifold  of  (reduced-order)  symmetric  positive 
definite  matrices  SPD(k)  [structure] 


ONLINE  INTERPOLATION  ON  A  MANIFOLD 


Given  a  p-parameter  system,  the  appropriate  Riemannian 
manifold,  and  its  logarithmic  and  exponential  maps 


SHOWCASE  APPLICATION 


F-16  Block  40 


CFD  model  FEM  structural  model 


Assisting  flight  test 

-  single  aircraft  configuration  - »  single  structural  ROM 

-  hundreds  of  flight  conditions  - -  database  of  fluid  ROBs 


FLIGHT  ENVELOPPE 


Mach  number 


ROM  DATABASE 


83  pairs  of  flight  conditions  (operating  points) 

-  70.6  hours  CPU  on  a  64-processor  Linux  cluster 


ON-DEMAND  PREDICTIONS 


Fast  responses  to  5  queries 

-  possible  scenarios:  flight  test,  flutter  clearance,  optimization 

-  interpolation  of  fluid  ROMs 


Lift  (Ibf) 


DELIVERED  ACCURACY 


Deep  transonic  point  (M^,  a)  =  (0.930, 1.3°) 


—  L-HFM 


Time  (sec) 


o 


1 


TURNAROUND  TIME 


/A 


F-16  Block  40  —  1  operation  point  —  1-processor  (desktop) 


CPU  time  for  interpolating  and  processing 
1  CFD-based  aeroelastic  ROM 


Fluid  ROM  interpolation  (5  points) 

0.2  second 

Aeroelastic  ROM  processing  (FD) 

0.3  second 

Total  CPU  time 

0.5  second 

>  real-time,  parameterized,  CFD-based  flutter  analysis 


IMPLEMENTATION 


MOBILE  DEVICES 


IMPLEMENTATI' 


CFD  Based  Acroclasli 


ON  MOBILE  DEVICES 


VISION 
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